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Recent advances in scanning tunneling spectroscopy performed on heavy-fermion metals provide 
a window onto local electronic properties of composite heavy-electron quasiparticles. Here we theo- 
retically investigate the energy and temperature evolution of single-particle spectra and their quasi- 
particle interference caused by point-like impurities in the framework of a periodic Anderson model. 
By numerically solving dynamical-mean- field-theory equations, we are able to access all tempera- 
tures and to capture the crossover from weakly interacting c and / electrons to fully coherent heavy 
quasiparticles. Remarkably, this crossover occurs in a dynamical fashion at an energy- dependent 
crossover temperature. We study in detail the associated Fermi-surface reconstruction and charac- 
terize the incoherent regime near the Kondo temperature. Finally, we link our results to current 
heavy-fermion experiments. 



I. INTRODUCTION 

Heavy-fermion metals where strongly localized / 
electrons tend to form magnetic moments which interact 
with more delocalized conduction (c) electrons, consti- 
tute an active and fascinating research area in condensed 
matter physics: Many non-trivial phenomena like com- 
peting orders, quantum criticality, unconventional super- 
conductivity, and quantum Griffiths phases find their re- 
alization in members of the heavy-fermion family.!^ De- 
spite significant theoretical progress in the field of cor- 
related electrons, the rich physics of heavy-fermion ma- 
terials remains only partially understood. The difficul- 
ties lie with spatially non-local phenomena - like exotic 
magnetism and superconductivity - while we have at 
least a qualitative understanding of the local process of 
temperature-dependent heavy-quasiparticle formation. 

Experimentally, most investigations have concentrated 
on thermodynamic or transport properties, which are ac- 
cessible in a straightforward fashion. In contrast, exper- 
imental results using single-electron spectroscopy in the 
heavy-fermion regime are scarce, as the required energy 
resolution in the sub-meV range is difficult to reach with 
present-day photoemission techniques. 

Spectroscopic-imaging scanning tunneling microscopy 
(SI-STM) is a surface-sensitive probe which allows to 
measure single-particle spectra with the required en- 
ergy resolution. Recent efforts in preparing clean 
surfaces of layered heavy-fermion compounds have al- 
lowed for the first time to measure spatially resolved 
maps of the differential conductance, dl/dV^ in heavy- 
fermion compounds .^"^ Thanks to impurity scattering 
processes, such imaging also allows to partially re- 
construct momentum-space information on the elec- 
tronic spectra, via energy-dependent Friedel oscillations, 
dubbed "quasiparticle interference" (QPI).^ Thus, SI- 
STM provides a unique opportunity to visualize the 
many-body quantum physics generated by the coupling 



between / and c electrons. The experimental data ob- 
tained on URu2Si^^ show a periodic lattice of Fano- 
shaped tunneling spectra at lowest temperatures, which 
is likely to arise from a combination of the hybridization 
between c and / bands and interference effects of differ- 
ent tunneling paths. In addition, the QPI signal shows 
signatures of heavy-band formation near the Fermi level. 

The theoretical modeling of these tunneling spectra 
has so far been restricted to a few easily accessible lim- 
its. Slave-boson mean- field techniques have been applied 
to the Kondo lattice model, which describe elements 
of the low-energy and low-temperature physics of heavy 
quasiparticles: The mean- field Hamiltonian is that of two 
hybridized bands of otherwise non-interacting fermions. 
As all inelastic scattering processes are neglected, phys- 
ical properties at elevated energies or temperatures can- 
not be described: Among the artifacts are a hard hy- 
bridization gap in the heavy-fermion state and an artifi- 
cial finite-temperature transition between this and a de- 
coupled high-temperature state. A model of heavy quasi- 
particle bands has been supplied by phenomenological 
Fermi-liquid self-energies in Ref. 16 in order to capture 
the partial filling of the hybridization gap. Such a quasi- 
particle broadening has been shown to be essential to the 
observation of the Fano-shaped peak in a Kondo lattice 
system. However, a consistent modeling for all energies 
and temperatures is not yet available. 

The purpose of this paper is to close this gap: We 
provide a detailed study of temperature-dependent elec- 
tronic spectra and QPI phenomena in a generic model 
of heavy-fermion metals, the periodic Anderson model 
(PAM). To this end, we numerically solve dynamical 
mean-field theory (DMFT) equations using Wilson's nu- 
merical renormalization group (NRG), a method which 
provides real-frequency spectra at all temperatures. The 
calculations allow us to track the formation of coherent 
heavy quasiparticles as function of energy and temper- 
ature, thus providing a basis for the interpretation of 
SI-STM experiments. As detailed below, we find that 
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the crossover temperature associated with the Kondo- 
driven Fermi-surface reconstruction is energy-dependent, 
demonstrating the dynamical character of the screening 
process. 



A. Outline 

The paper is organized as follows. In Sec.|TT]we intro- 
duce the PAM and its treatment within DMFT using the 
NRG impurity solver. Sec. |III| describes the calculational 
scheme for the tunneling and QPI signals, the latter in- 
volving the scattering off isolated impurities. In Sec. IV 
we discuss various general aspects of heavy fermions and 
their description within a local self-energy approxima- 
tion, with focus on energy scales and corresponding fea- 
tures in the single-particle spectra. Sec. |V| is devoted 
to a detailed discussion of our numerical results, where 
we present both momentum-integrated and momentum- 
resolved single-particle spectra as well as constant-energy 
cuts through single-particle and QPI spectra, all as func- 
tion of temperature. This will allow for a detailed dis- 
cussion of the crossover from weakly interacting c and / 
bands at high T, with a small Fermi surface, to coherent 
heavy quasiparticles with a large Fermi surface at low T. 
A summary of experimentally relevant aspects closes the 
paper. 

Readers mainly interested in the temperature-driven 
band reconstruction from light to heavy quasiparticles 
should jump ahead to Sec. |VB| and Fig. [6j 



II. MODEL AND DYNAMICAL MEAN-FIELD 
APPROXIMATION 

Here we describe the microscopic approach taken to 
calculate single-particle spectra. The concrete modeling 
of SI-STM data is covered in Sec. IIIII 



is the hybridization between the two fermion bands, and 
U is the local Coulomb repulsion between the /-electrons, 
which is often to be taken the largest energy scale. 

For negative e/ and large [/, local moment tends to 
form on the / orbitals. The charge fluctuation scale of the 
/-electrons is determined by the Anderson width Fq = 
(l/-^)Ek^k^(^k). Provided that Fq <C /7, a mixed- 
valence regime is reached for |e/| ~ Fq with < n/ < 1, 
whereas |e/| Fq leads to stable local moments. In this 
regime, charge fluctuations can be integrated out, and 
one obtains the Kondo lattice model 



(2) 



kcr 



where the impurity spin Si is coupled to the conduction 
electrons spin density at site i, Si = ^^^f ^i^jTcja'Cia' l'^ 
and = (l/A/") Xlk ^^'^^^^kcr- For a local hybridization, 
Vk = the Kondo coupling Jk in Q is related to the 
parameters of the Anderson model through 



Jk = 



1 



1 



(3) 



The Anderson model ([T]) describes orbitally non- 
degenerate / states - this situation often applies to Ce- 
based heavy-fermion systems: the lowest / conflguration 
is /^, whose multiplicity is reduced to a Kramers dou- 
blet due to crystalline-electric-fleld (CEF) splitting, and 
which is well separated from the /^ state. The situation 
is more complicated, e.g., in uranium compounds where 
/^ conflgurations cannot be neglected. It is, however, 
believed that Eq. ([T]) captures, at least qualitatively, the 
low-energy physics of many heavy-fermion compounds. 



B. Dynamical mean- field theory 



A. Periodic Anderson model 

A standard microscopic model used in the context of 
heavy-fermion compounds is the so-called periodic An- 
derson (or Anderson lattice) model, ^ which describes hy- 
bridized c and / bands with a strong local repulsion on 
the / orbitals: 

^ = E hcLck. + e//L/k<. + Vk(cL/k. + H.c.)" 

kcr 

Here cj^cr(/kcr) create a conduction electron (/-electron) 
of momentum k and spin a, and = fj^fia counts the 
/ particles at site i. ek and e/ are the band energies of the 
conduction and dispersionless /-electrons, respectively, 
measured relative to the chemical potential. Finally, Vk 



As the periodic Anderson model with non-zero U is 
not exactly solvable, further approximations have to be 
made. A well-established and successful method, which 
is able to capture the effects of strong local correlations, 
is the so-called dynamical mean-fleld theory.'^^'^^ Here, 
the self-energy due to the Hubbard-like interaction U is 
assumed to be momentum-independent, 

S/(k,c^)^SyH. (4) 

This approximation becomes exact in the limit of infl- 
nite coordination number. As a result of Eq. Q, the 
lattice problem Eq. ([T]) can be mapped onto an effective 
single-impurity Anderson model supplemented by a self- 
consistency condition for the effective medium (or bath) 
of the impurity. Assuming that both c and / electrons 
live on the same lattice structure, with a local hybridiza- 
tion Vic = the self-consistency equation for the single- 
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particle propagator reads 



M ^ iuj - Ef - S/(iw) 



V2 



(5) 



de- 



Po(e) 



1 

iuj — Cf — A{iuj) — Tif{iuj) 



y2 

iu; — e 



^SIAM 



where Q^j{iuj) is the local / Green's function, with 
the superscript denoting the impurity-free system (see 
below), and po(e) = (l/A/") Xlk ^(^ ~ ^k) is the non- 
interacting c electron density of states (DOS). Further, 
A{iuj) is a generalized hybridization function, which de- 
fine the impurity problem and depends implicitly on the 
/ self-energy T,f{iuo) and is thus different from its non- 
interacting form due to the / correlations effects. 

Using the / self-energy one can express the 
momentum-dependent Green's functions as follows: 



, (6) 
, (7) 



{iuj — e\^){iuj — ef — T^f{iuj)) — 
iuj — Ck 

{iuj — e\^){iuj — — Tif{iuj)) — ' 
y 



Here, ^^^(k, icj) is the Fourier transform of ^^^(k, r) = 

— (T7-Ckcr('?")C][^cr)' with Tt- being the time-ordering opera- 
tor on the imaginary time axis; Sjj(k, iuo) and Q^fi}^-, ioj) 
are defined analogously. Note that all Green's functions 
are independent of spin in the considered paramagnetic 
state. From Eq. (|6|, we can define a self-energy Sc(ico') 
for the conduction electrons according to 



Sc(^Cj) 



IUJ ■ 



(9) 



C. Numerical renormalization group 



[±V^A-(^+i) , ±V^A-^] forn = 0, 1, . . . , oo, where W is 
the half-bandwidth of the bare conduction band. The 
original Hamiltonian is then mapped onto a semi-infinite 
chain, the Wilson chain, where the Nth link represents 
an exponentially decreasing energy scale ^ A~^/^. The 
chain Hamiltonian is diagonalized iteratively, by starting 
from the impurity and successively adding chain sites. In 
each step, the high-energy states are truncated to main- 
tain a manageable number of states Ns. The reduced 
states are expected to capture the spectrum of the origi- 
nal Hamiltonian on a scale ~ A~^/^, corresponding to a 
temperature T = W/ PA~^^'^ from which all thermody- 
namic expectation values are calculated. Spectral infor- 
mation at the temperature T is calculated by collecting 
information from NRG steps 1, . . . , A/" which yields dis- 
crete spectra on a logarithmic energy scale down to ener- 
gies of (a fraction of) T. The impurity self-energy needed 
in the DMFT iteration is most accurately evaluated as 
the ratio of two propagators, eliminating the need for us- 
ing the Dyson equation. For more details of NRG and 
its application to DMFT we refer the reader to a recent 
review article.'!^ 

The DMFT-NRG method has been employed to inves- 
tigate heavy-fermion physics in the past.^^ 26 29 jJq^_ 
ever, those investigations mainly focused on thermody- 
namic properties or momentum-integrated spectra, and 
consequently the Bethe lattice was used. Here, we are 
motivated by layered heavy-fermion materials being good 
candidates for STM investigations, and thus we per- 
form all numerical calculations for a two-dimensional 
(2d) square lattice with nearest-neighbor hopping. NRG 
parameters are A = 2, Ng = 600, and P = 1 unless oth- 
erwise noted. 



III. TUNNELING SPECTRA AND 
QUASIPARTICLE INTERFERENCE 

After having described our theoretical approach to the 
clean (i.e. impurity-free) bulk heavy-fermion system, we 
now turn to aspects relevant to SI- STM. 



To solve the effective impurity problem arising within 
DMFT or its generalizations, a variety of different meth- 
ods have been developed, all with individual advantages 
and drawbacks. In the present case, a non-perturbative 
method is preferable which (i) can access the small en- 
ergies and temperatures relevant for heavy-fermion for- 
mation and (ii) directly provides spectral function on the 
real frequency axis, such that analytic continuation from 
the imaginary axis is not required. 

Here, we choose Wilson's numerical renormalization 
group (NRG) technique, which has been successfully 
implemented in the context of DMFT for the Hubbard 
model, for periodic Anderson model, and the Kondo 
lattice model.'^^''^ NRG is based on a logarithmic dis- 
cretization of the energy axis, controlled by a param- 
eter A > 1: The energy axis is divided into intervals 



A. Tunneling spectra 

In the context of tunneling experiments both on sin- 
gle magnetic impuritie s on m etallic surfaces^^ and on 
heavy-fermion systems,'^^'!^ it has been proposed that 
the tunneling current arises from the interference of two 
channels, namely a tunneling process from the tip - po- 
sitioned over a particular site - into a conduction- 
electron state, with an amplitude tc, and another one 
to an /-electron state, with an amplitude tf. Assuming 
both processes to be spatially local, the tunneling piece 
of the Hamiltonian may be written as: 

Ht = J2 {^ccl + tffl) p. + H.C., (10) 
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where Pa- destroys an electron with spin a in the tip. 

Assuming further that the tip and the system are in 
thermal equilibriump^ the total tunneling current, flow- 
ing from the tip into the system, to lowest order in the 
tunneling amplitudes t/,tc, is given by 

2e 



/(r„ eV) = ^ J du ptip(c^ - eV) [f{u - eV) - f{u)] x 



(11) 



Here, ptip is the tip DOS, V the applied bias voltage, 
f{oj) = 1/(1 + e^^^) the Fermi-Dirac function, and all 
real-frequency Green's functions are to be understood as 
retarded ones. In a translation-invariant system, the real- 
space Green's function Gcd'^i^ ^j^^) depends on (r^ — r^) 
only, such that /(r^, eV) is independent of the tip position 
r^. (Note that we do not account for the sub-atomic 
structure of the Bloch states.) The last term in (11) 



describes the quantum interference processes between the 
two tunneling paths. 

If the tip DOS Ptip(^) is independent of a;, the differ- 
ential tunneling conductance reduces to 



dl 2e 



du f'{u - eV) X 

+ 2tet/ac/(r^,r„cj)],(12) 
which is the quantity of interest in discussing STM data. 



B. Impurities 

Friedel oscillations and QPI phenomena rely on the 
breaking of translational invariance due to impurities or 
crystal imperfections which act as scattering centers. In 
fact, intentional impurity doping has been employed to 
enhance QPI signatures in SI-STM experiments.^ 

As we are interested in qualitative features, we use the 
simplest model capable of describing interference phe- 
nomena in the framework of the Anderson model. Con- 
sider a random distribution of Ug non-magnetic point-like 
scatterers at sites R^, z = 1, . . . , ng which couple to the 
conduction electrons only. This introduces a scattering 
term Hs- in the Hamiltonian 



Us = QsY. 4(R^)Ca(R^ 



= l,cr 



i=l,cr 



(13) 



where Qs is the strength of the impurity potential and 
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b) 
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FIG. 1: a) Feynman diagrams which enter the single-impurity 
T-matrix used in our calculation. Bold lines are full (i.e. inter- 
acting) propagators the cross denotes impurity scattering 
with strength Qs. b) Sample diagrams contained in a), where 
thin lines are now non-interacting propagators, and the wavy 
line represents the interaction U . c) An interference process 
between Qs and U which is neglected in a). 



Apparently, more general forms of Qs are possible, but 
will not be considered here. 

For static point-like impurities embedded into a free- 
electron gas {U = in our case), the Green's function 
matrix Q in the presence of scattering can be expressed 
exactly via the T-matrix as follows: 



g{vy,uj) = g°{r-r',uj) 



(15) 



with 



1=1 



and the single-impurity T-matrix i{uj) 



(16) 



(17) 



where 1 is the identity matrix in spatial, spin and fermion 
indices, and now denotes Green's functions in the ab- 
sence of scattering. 

In our calculations for the interacting system, we shall 
employ the following approximations: (i) We neglect in- 
terference process between different impurities, formally 
Tij idjj . (ii) We capture the scattering effects using 
Eqs. (15) and ([Tt]), i.e., we evaluate a single-impurity 



T-matrix t according to (17) with full (i.e. interacting) 



propagators Q^. This neglects processes where impurity 
scattering and electron-electron interactions interfere, as 
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illustrated in Fig.[T] The change in local density of states 
induced by the potential scattering can then be expressed 
as 



1 



(18) 



-ImTr 



(iii) Unless otherwise noted, we employ the lowest-order 
Born approximation, i.e., t = Qs- 

The approximations (i)-(iii) are designed to capture 
the effect of weak impurity scattering. They do not ac- 
count for impurity- induced local changes in the Kondo 
screening j ^^ l ^^ l a description of which would require a 
solution of real-space DMFT equations with a self- 
consistent set of impurity problems^^ - a task which we 
leave for future work. 



C. Quasiparticle interference 



If we Fourier transform (5p(r, uj) ( 18 ) into momentum 



space and restrict ourselves to the lowest Born approxi- 
mation, we find 



(5p(q,cj) 



-1 



E 

k 



^0(k,w)5,^0(k + q, 



(19) 

i.e, to leading order in the impurity strength, the FT of 
the LDOS separates into a factor from the QPI of de- 
scribing "band structure", and another factor represent- 
ing the spatial distribution of the scatters. The latter is 
a smooth function of momentum; it is relevant for the 
overall intensity distribution of 5p(q, cj) in momentum 
space, but not for the presence of sharp features (peaks, 
ripples etc.). 

Therefore, we may restrict the calculation to the case of 
a single impurity, say at Ri = 0. We note that this type 
of approximation (single impurity, T-matrix approxima- 
tion in the presence of interactions) has been frequently 
used to describe QPI phenomena in various physical sys- 
tems in the past (see e.g. Refs. I12|13|36|37|) . Using 
(12), the impurity-induced change in differential tunnel- 
ing conductance, i.e. its spatially inhomogeneous piece, 
is given 



2e 

-QsPtipj^ I duj f'{uj - eV) X 



im [tlgt,{v,Lo)gt,{-v,uj) + t}g%{v,uj)g%{-v,uj) 

+t,tfglir,u;)g°fi-r,u,) + t,tfg%{v,Lo)gl{-v,^) 



(20) 



Note that no appears here because of the structure of 
the scattering potential Qs (14). Transforming to Fourier 



space and assuming inversion symmetry, we obtain 
5^(q, eV) cx ^ Im [tl g°,(k, el^)g°,(k + q, eV) 

k 

+ t2 gO^(k^eF)gO^(k + q,eT/) 

+ 2tet/ a°e(k,ey)a°/(k + q,el^)] 

(21) 

which is the equation we will use below to generate nu- 
merical results describing QPI. Note that the Fourier- 
transformed LDOS is complex in general, but the present 
situation obeys inversion symmetry w.r.t. the position of 
the single impurity, and hence dI/dV{q_^eV) is real. 

In general, QPI phenomena as captured, e.g., by equa- 
tions (18)-([2T]), lead to energy-dependent Friedel oscil- 



lations in the LDOS, caused by scattering off the im- 
purities. Those oscillations at an energy u are primar- 
ily determined by the shape of the iso-energy surface at 
UJ of the underlying band structure, i.e., the oscillation 
wavevectors q are given by approximate nesting wavevec- 
tors or by wavevectors which connect points of high DOS 
in momentum space. (Note that this argument neglects 
the infiuence of the real parts of the propagators which 
enter in Eqs. (19)-(21) as well.) Thus, the Fourier trans- 



form of the LDOS "ripples" observed in SLSTM ex- 
periments allows to approximately extract characteristic 
wavevectors q(cj) of the underlying electronic state. Be- 
low we shall show that QPI in heavy-fermion systems 
is capable of detecting the Fermi-surface reconstruction 
from high to low temperatures. 



IV. ENERGY SCALES AND GENERAL 
CONSIDERATIONS 

Before diving into numerical results, let us discuss a 
few general aspects of heavy-fermion physics relevant to 
SI-STM experiments performed at variable temperature. 

Provided that the Anderson model ([T]) is close to the 
Kondo limit, i.e., /7, — e/ > V'^ /W where W is the half- 
bandwidth of the c electrons, local moments will form 
on the / sites upon cooling below a temperature of order 
U. At those high temperatures, conduction electrons will 
weakly scatter off those local moments, a process which 
can be described in perturbation theory in the Kondo 
coupling Jk in ([2|. Upon cooling, the scattering intensity 
becomes large when T reaches the single-impurity Kondo 
temperature Tk where perturbation theory breaks down. 

In the opposite limit of low temperatures (and in the 
absence of spontaneous symmetry breaking e.g. by mag- 
netism or superconductivity) , a heavy Fermi liquid forms 
away from half-filling. This Fermi liquid is adiabatically 
connected to the weakly interacting limit of the PAM 
([T]), and perturbation theory in U is at least formally 
applicable. The Fermi liquid is bounded in temperature 
from above by a coherence temperature Tcoh which acts 
as an effective Fermi energy for the heavy quasiparticles. 
Typically Tcoh < Tk- 
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slave boson 6^ is introduced such that 



2 

ky — kv 



FIG. 2: Schematic representation of the T ^ band struc- 
ture expected for a PAM. Feff and e/ are the renormalized 
hybridization and f-chemical potential respectively. is the 
hybridization gap between the lower and upper renormalized 
bands. 



In the low-temperature Fermi-liquid regime, Lut- 
tinger's theorem applies, and consequently the Fermi 
volume is "large", i.e., given by Vfl = Kd{ntotmod2)^ 
where ntot = + n j is the total number of electrons per 
unit cell and = (27r)^/ {2vo) a phase space factor, with 
^0 the unit cell volume and d the spatial dimensionality. 
In the limit of small Tk, it also makes sense to consider 
a Fermi volume for temperatures Tk <C T <C W^. Here, 
the / electrons do not contribute, leading to a "small" 
Fermi volume with Vg = KdUc. 

Inelastic scattering between c and / electrons is 
particularly strong for intermediate temperatures near 
^coh,^K- In fact, a common experimental definition of 
Tcoh is via the maximum in the resistivity p{T). Thus, 
one expects that the temperature-driven crossover be- 
tween the small and large Fermi volumes occurs at a 
temperature of order Tcoh,^K- 

The physics of the heavy Fermi liquid can be under- 
stood, to some degree, in a picture of renormalized bands: 
Two bands of non-interacting c and / fermions hybridize 
with a renormalized hybridization Vq^. This results in 
two quasiparticle bands with dispersions 



(23) 



± _ ek + e/ ± V(ek - £/) 



■ 4v;l 



(22) 



In the simplest approximation, h is uniform and static 
and measures the effective hybridization between the 
bands, i.e., Feff in Eq. (22) becomes ]4ff = (Kff = 
6Jk) in the Anderson (Kondo) model. (Note that h 
is also a measure of the mass renormalization of the 
quasiparticles, i.e., the mass enhancement is given by 
m*/m = 1/^^.) This mean- field approximation can be 
formally justified in a limit where the spin symmetry of 
the original model is extended from SU(2) to SU(A^) with 
^ oo. Remarkably, the qualitative picture of effective 
quasiparticle bands. Fig. [2j has been reproduced in nu- 
merical treatments of the periodic Anderson (or Kondo 
lattice) models, using both DMFT^^^^ and its cluster 
generalizations^^. Hence, a quasiparticle description of 
low-temperature heavy- fermion bands appears justified, 
with a number of caveats noted in the following. 

While the slave-boson mean-field theory correctly cap- 
tures the exponential dependence of the Kondo scale Tk 
on Jk, W ^ and accounts for the enlargement of the Fermi 
surface in the Fermi-liquid regime, it misses all inelas- 
tic scattering processes, as it only operates with non- 
interacting fermions. Consequently, the physics not only 
at finite temperatures, but also that at T ^ and finite 
energies is not described well. Both interaction-induced 
broadening of spectral features away from the Fermi level 
and incoherent scattering at elevated temperatures are 
not included. For instance, the hybridization gap in the 
heavy-fermion state. Fig. [2j is predicted to be a sharp 
gap in the mean-field theory, but can be expected to 
be significantly smeared or even washed out by inter- 
action effects. In fact, it has been argued that the exper- 
imentally observed Fano-line shape is a consequence of a 
quasip articl e broadening not captured by the mean-field 
theory.^^I More seriously, the crossover from T <C Tk 
to T Tk turns into an artificial phase transition at 
Tk in the slave-boson theory. Hence, to understand the 
crossover between large and small Fermi volumes upon 
varying temperature requires at least to treat local corre- 
lation effects in a non-perturbative manner - this is what 
we shall do below using DMFT. 



where e/ is the renormalized energy of the /-electrons. 
These bands describe sharp quasiparticles formed as a 
mixture of / and c degrees of freedom. The renormalized 
bands, shown schematically in Fig. |2j are separated in 
energy by a direct (or optical) gap Aopt = 2V^ff at the 
crossing of the original conduction band and e/, while 
the bottom of the upper band is separated from the top 
of the lower band by an indirect gap (x \/V^^W . Ag 
is typically called hybridization gap. 

While this two-band picture is trivially realized at = 
in ([T]), with V^ff = V and if = e/, it can be obtained 
using a slave-boson mean-field approximation from either 
Eq. ([T]) in the U ^ oo limit or from Eq. ([2|. In both 



V. NUMERICAL RESULTS 

We have performed extensive studies of the paramag- 
netic phase of the PAM using DMFT-NRG at various 
temperatures, tracking the formation of the heavy Fermi 
liquid upon cooling out of the light conduction band and 
the /-electrons. 

For simplicity, we consider a square lattice with 
nearest-neighbor hopping. 



ek = — 2t(cos kx + cos ky) + Cc 



(24) 



and half-bandwidth W = 4t which we employ as energy 
unit. We assume the chemical potential to be located at 
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FIG. 3: Frequency dependence of the self-energies, ImS/ and 
ImEc, of the / and c electrons, respectively, calculated using 
DMFT-NRG for the periodic Anderson model with parameter 
U — 10,6/ — — 0.9, Cc — 0.4, y — 0.55 at a low temperature 
of T = 1.8 X 10~^. Energies are given in units of the half- 
bandwidth W. 



zero energy, such that Cc controls the conduction-band 
filling. 

Most of the figures shown have been obtained with the 
following set of the PAM parameters: Cc = 0.4, = 
-0.9, V = 0.55, U = 10 (in units of W). With these 
parameters the Anderson width is Fq ~ 0.51. In the 
T ^ limit, we find the occupation number of the c (/) 
electrons is Uc ~ 0.57 {uf ^ 0.87), i.e., we work close to 
the Kondo regime, with a deviation from integer filling 
comparable to that in actual Ce or Yb heavy-fermion 
systems. The temperature variation of the band fillings 
is less than 5% within the temperature range considered 
here. For the above parameters, we find a T ^ mass 
renormalization 



m*/m 



duj 



60 



(25) 



which allows to define a low-temperature Kondo scale 
7^ = Fq m/m* ^ 8 x 10~^ - this roughly matches the 
full width at half maximum (FWHM) of the T = 
Abrikosov-Suhl resonance in see Fig. [5] below. Note, 
however, that a unique definition of neither the Kondo 
nor the coherence scale exists, and depending on ones 
choice may differ by up to an order of magnitude. In 
our case, we find the temperature scale where the impu- 
rity entropy within DMFT reaches 0.5 In 2 to be T1/2 ^ 
1.5 X 10~^, and heavy-fermion bands are fully formed 
only below a coherence temperature of Tcoh ~ 2 x 10~^. 

For the 2d system considered here, the van- Hove singu- 
larity (vHs) arising from the saddle point in the disper- 
sion is a prominent feature in the DOS which will play a 
role in all momentum-integrated spectra discussed below. 



A. Local electronic spectra and self-energies 

Once we have solved the self-consistent equation ([5| for 
the / self-energy Sj, we can form the Green's function 
matrix whose elements are given in (pHsl). In this sec- 
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FIG. 4: a) c and b) / electron local spectral functions of the 
periodic Anderson model, calculated using DMFT-NRG at 
T = 1.8 x 10-^ The different curves correspond to different 
values of the local interaction [/; other parameters are as in 
Fig. [3] Because the chemical potential is kept fixed, the total 
and individual band fillings vary as follows: ntot — ric — rif ~ 
2.2 - 0.42 - 1.78 {U = 0), 1.98 - 0.78 - 1.2 {U = 1.5), 1.52 - 
0.63 - 0.89 {U = 3.5), and 1.44 - 0.57 - 0.87 {U = 10). 



tion we discuss self-energies as well as the local spectral 
functions the c and / electrons at different temperatures. 



1. Self- energies 

As a check of Fermi-liquid behavior, we start with the 
frequency dependence of the imaginary parts of the low- 
temperature electronic self-energies, ImSc and ImSj, 
displayed in Fig.jSj Both exhibit quadratic behavior in cj, 
in a limited energy range around the Fermi level, as ex- 
pected for a Fermi liquid. We note that the self-energies 
have a small residual value at = 0, arising from both in- 
accuracies of NRG in fulfilling spectral sum rule^^ and 
the artificial broadening which has been introduced to 
stabilize the numerics. (Our results do not depend on 
details of this broadening.) 



2. Low-temperature spectra 

Next we turn to the local spectral functions, Pc(^) = 
( — l/7r)Im^cc(^) and pfioo) = {—l/7T)lmQff{uj). To il- 
lustrate the effect of interactions and the concomitant 
deviations from the two-band picture of slave-boson the- 



ory explained in Sec. [TV] we show the U dependence of 
the low-temperature spectra in Fig. |4j For the c spectral 
function. Fig. |4k, increasing U has the following effects: 
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FIG. 5: Local spectral functions a) pc{(^) and b) pf{uj) of the PAM calculated using DMFT-NRG for different temperatures 
T, with model parameters as in Fig. [s] al,bl) Spectra on the full energy range for two temperatures, T ^ Tk and T <^ Tk- 
a2,b2) Temperature evolution of pc and pf shown in a restricted energy range near the Fermi level. a3,c3) Low-temperature, 
low-energy results for pc and p/, with the high- temperature spectra at T — 4.7 x 10~^ subtracted. The progressive formation 
of the hybridization pseudogap (panel a) and the Abrikosov-Suhl resonance (panel b) are clearly visible. 



a renormalization of the band positions arising from the 
real part of the self-energy (note that the chemical po- 
tential is kept fixed), a renormalization and simultane- 
ous displacement of the inter-band gap, and its smearing 
due to the quasiparticle broadening induced by ImSc(co'). 
While the two first effects are accounted for in the slave- 
boson theory, the last is absent in this approximation. 
The progressive smearing of the hard gap at /7 = and 
its displacement are also visible in the / spectrum. Fig. 
^p. Increasing U induces the formation of the Kondo 
peak (or Abrikosov-Suhl resonance) near the Fermi level. 
Its width can be considered as a measure of the Kondo 
scale Tk, which approaches a finite value in the U ^ oo 
limit (due to the finite e/). 

A few comments are in order: A sharply defined hy- 
bridization gap in corresponding to the indirect gap 
of Fig. pj is only present for U ^2, while it is re- 
duced to a dip or pseudogap for larger values of U ^3. 
Remarkably, with increasing U one observes an appar- 
ent pinning of the vHs of the c band to the Fermi level. 
As will become clear from the band structures discussed 
below, this is due to the progressive narrowing of the 
lower quasiparticle band at low energies. As a result, the 
vHs and hybridization gap feature cannot be clearly sep- 
arated in pc (both exist on an energy scale of order Tk). 
The vHs is also responsible for the strong particle-hole 
asymmetry of the c spectra. In pf both the vHs and the 
hybridization gap are less prominent. Finally, we remind 
the reader that our single-orbital PAM does not account 
for multiple / levels and their crystal- field splitting, i.e.. 



we only model the physics of the lowest-energy crystal- 
field doublet. 



3. Temperature evolution of spectra 

The temperature evolution of both c and / local spec- 
tral functions in the Kondo regime {U = 10) is shown 
in Fig. |5] As can be seen in Fig. ^jpl, the / spectrum 
has the well-known three-peak structure at low temper- 
ature: the two charge excitation peaks at e/ and e/ + /7, 
well separated from the Abrikosov-Suhl resonance at the 
Fermi level. The latter is absent at high temperature 
and is the fingerprint of the Kondo effect. A close-up 
view of Pf{uj) near the Fermi level. Figs. ^2 and |5]33, 
shows the development of this resonance with decreasing 
temperatures. 

In the local c spectrum. Fig. [5^1, cooling induces the 
hybridization pseudogap and a second van-Hove singu- 
larity near the Fermi level. This is again shown in de- 
tail in Figs. [5^2 and [5^3, where the development of a 
positive-frequency dip with decreasing temperatures is 
emphasized. As noted above, this dip is what remains 
of the hybridization gap once inelastic scattering is fully 
taken into account; its particle-hole asymmetry is a band- 
structure effect arising from the low-energy vHs in the 
present 2d case. The displacement of the higher-energy 
vHs in Fig. [5^2 can be taken as a measure of ReSc. The 
low-energy structures apparent at intermediate temper- 
ature ^ 10"^ . . . 10"^ are connected with the dynami- 
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FIG. 6: Evolution of the electronic spectrum of the PAM, with parameters as in Fig. |3] with temperature: Shown is 
(— l/7r)Im(^// + ^cc) (kjCj) as function of kx = ky and u over a) the full range of energies and b) a restricted energy win- 
dow near the Fermi level. In a) the dotted horizontal lines mark the renormalized / level (i.e. the lower Hubbard band), 
while in b) the dotted lines represent the Fermi level. Further the dashed line in al) shows the bare conduction band ek- At 
low temperatures, the formation of heavy bands is obvious, with the reconstruction happening mainly within the optical gap. 
In a5) the dashed lines represent a fit of the intensity at elevated energies to the effective two-band picture using (22) with 
Veii — 0.17, €f — 0.025, and ec = 0.47, for detail see text. 



cal character of the band reconstruction and will be dis- 
cussed in more detail below. 



B. Renormalized band structure 

After presenting the local (i.e. momentum-integrated) 
spectral functions, we turn to the full momentum de- 
pendence of the spectra. (Recall that in DMFT 
the interaction-induced self-energies are momentum- 
independent, but the Green's functions acquire mo- 
mentum dependence from the non-interacting bands.) 
Fig. [6] illustrates the temperature evolution of the elec- 
tronic spectrum in a false-color plot of (— l/7r)Im(^jj + 
^cc)(k, uj) as function of energy uo and momentum k along 
the line kx = ky. At high temperature T ^ T^^ Fig.|6^1, 
the signal consists of the rather sharp bare c band and 
the "atomic" / levels at ej and Cf -\- the latter being 
broadened by the hybridization. In the opposite limit 
T <C Tx, Fig. [6^5, we have two self-energy-broadened 
bands, both being rather flat near the Fermi level. One 
may fit this low-temperature band structure using the 



two-band picture described in Sec. IV However, we find 
that there is no unique fit for all energies: The features at 
elevated energies are well described by the fit shown in 
Fig. [6^5, which, however, overestimates the low-energy 
slope of the heavy band by a factor of 2.5. Alternatively, 
a low-energy fit results in V^f^ = 0.072, consistent with 
/m ~ 60 obtained from DMFT self-energy - this im- 
plies that the DMFT result has to be understood in terms 
of an energy-dependent band hybridization. 



The reconstruction of the band structure from a weakly 
interacting c band and / levels at high temperature 
to two heavy bands at low temperature occurs mainly 
within the optical gap ~ 2]4ff, as can be seen from 
Fig. ^jp: Portions outside this window show little varia- 
tion with temperature. When temperature is decreased, 
spectral weight is transferred gradually from the bare c 
band near the small Fermi surface to the parts of the in- 
cipient heavy-fermion bands near the boundaries of the 
momentum window. Importantly, for a given energy this 
weight transfer happens at an energy-dependent temper- 
ature, i.e., upon cooling higher energies are reconstructed 
first. Figs. [6]32-[6]34. As a result of this dynamical re- 
construction, there is an "island" of spectral intensity 
near the bare c band left at intermediate temperature, 
Fig.|6]32, which shrinks and finally disappears upon cool- 
ing to significantly lower T. 

Constant-energy cuts of the spectral functions, now 
separated into (— l/7r)Im^jj and (— l/7r)Im^cc7 illustrat- 
ing the weight distribution in 2d momentum space, are 
shown in Fig. [t] for two different temperatures (higher 
than the one at which the renormalized bands are fully 
constructed, which is roug hly 2 X 10-"^). The plots al- 
low to track in detail the spectral-weight transfer from 
the high-temperature bare conduction band to the low- 
temperature renormalized bands; note that we have em- 
ployed a logarithmic color scale in order to visualized 
weak-intensity features as well. 

Starting the discussion with c electron spectrum, we 
confirm that there is essentially no temperature depen- 
dence for energies outside the optical gap, compare pan- 
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FIG. 7: Constant-energy cuts through the electronic spectra of the PAM, with parameters as in Fig.js] Panels a) and c) show 
(— l/7r)Im^cc(k, cj) in a quarter of the first Brillouin zone, panels b) and d) show ( — l/7r)Im^//(k, cj)? The temperatures are 
a,b) T — 7.3 X 10~^ and c,d) T = 1.8 x 10~^. Note that the intensity color scale is logarithmic in all cases. 



els al) and a6) to cl) and c6) in Fig. 7| In contrast, 
for energies inside the optical gap one observes a signif- 
icant temperature dependence: At higher temperatures 
a rather diffuse signal is present in panels a2)-a5), while 
quasiparticle peaks form at low T and negative energies, 
panels c2) and c3), and spectral weight disappears at low 
T and positives energies, panels c4) and c5) - the latter 
correspond to energies inside the hybridization gap. Note 
that panel Fig. [7^4 contains shows a well-defined (albeit 
weak) iso-energy contour - this represents a yet unrecon- 
structed piece of the light band which only disappears at 
lower temperature, see also Fig. (6)3. 

The / electron spectra in all panels. Fig. [7]3,d, show 
a pronounced weight reduction which follows closely the 
"small" Fermi surface defined by a; — ek — ReHcl^) = 0; 
this effect can be easily deduced from Eq. ([t]) (given that 
the interaction-induced self-energy has a non-zero imag- 
inary part). A well-defined quasiparticle peak is only 
visible in Qff at frequencies where portions of the new 
renormalized bands are constructed, i.e., in panels d2) 
and d3). 



C. Differential tunneling conductance 

In the presence of lattice translational symmetry, i.e., 
in the absence of impurities, the differential tunneling 
conductance dl / dVivi^uj)^ Eq. (12), is independent of 



the spatial position (i.e. lattice site) r^. It is a measure 
of the local density of states, but depends on the ratio 
tf /tc of the tunneling amplitudes. 

Fig.jsjshows our results for the energy (or bias- voltage) 
dependence of dl/dV for different values of tf/tc and 
each for three representative temperatures, T ^ Tk, 
T ~ Tx, and T <C Tk- For small values of tf/tc^ sig- 
natures of the hybridization gap are clearly seen, but - 
as expected - there is never a hard gap near the Fermi 
level, in contrast to results from slave-boson mean-field 
theory. ^^ ^^ The reason of course are inelastic scatter- 
ing processes, captured by DMFT but absent in the 
static mean-field approximation. (Those were added phe- 
nomenologically to explain the smearing of the hard gap 
near the Fermi level and the appearance of the asymmet- 
ric peak in experiments.^^) 
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FIG. 8: Evolution of the differential conductance dl/dV for 
small bias voltage with the ratio t//tc for different tempera- 
tures T > Tk, T - and T < Tk. The effect of destruc- 
tive [constructive] interference for positive [negative] ratios 
tf /tc is visible at low temperatures in panels b), d), f) [c) and 
e)], while the interference is essentially ineffective at high T. 



For = 0, Fig. [8^, the differential tunneling con- 
ductance is proportional to the c-electron LDOS. It is 
almost flat at high temperature while it has a broad 
peak and a dip at low T, resulting from the interplay 
of the vHs and the hybridization gap. The particle-hole 
asymmetry is thus inherited from the conduction band. 
This asymmetry appears opposite to the one observed 
experimentally^ in URu2Si2; in our calculations, such 
asymmetry would result upon choosing Uc > 1. (Note 
that the asymmetry in experiment may be influenced by 
many factors: conduction band dispersions and fillings, 
relevant crystal field levels of the Kondo ion, and the 
tunneling paths relevant for STM.) 

For tf /tc ^ 1, the main contribution to the differential 
tunneling conductance comes from the /-electron LDOS. 
Accordingly, a large (Kondo) peak is observed in dl/dV 
near a; = and the low-energy particle-hole asymmetry 
is inverted. This is shown in Fig. [SfD. 

For intermediate values of tj/tc, quantum interferences 
between the two tunneling channels have dramatic effects 
at low temperature as can be seen in Fig. [Sj^-f. These 
effects depend on the sign oitf /t^ With increasing pos- 
itive tf/tc and at low temperatures, the vHs peak at 
negative bias gradually decreases. Then the low-energy 



particle-hole asymmetry in dl/dV is inverted. Fig. [8]i,f, 
while the Kondo peak emerges near u = 0. Notice that, 
while the ratio tf/tc increases, the differential tunneling 
conductance at low T decreases below the one at high 
temperature before it increases drastically for t//tc ^ 1. 
The initial decrease is due to destructive interference 
of the two tunneling channels arising from a negative 
Im^c/(^) (for the uj range considered here) and positive 
ratio tf/tc- (This effect is essentially ineffective at high 
temperature.) Conversely, for negative tf/tc we observe 
constructive interference of the tunneling channels at low 
UJ and T, as can be seen in Fig. [8]3,e. Note that the type 
of interference (constructive vs. destructive) depends on 
the sign of tftclraQcfiuo)^ the latter in turn depending on 
particle-hole asymmetry and microscopic details. 



D. Quasiparticle interference 

Finally, we study quasiparticle-interference phenom- 
ena induced by scattering off impurities. As detailed in 



Sec. |IIIB] we concentrate on the physics of a single point- 
like scatterer in the Born limit, which should be appropri- 
ate for dilute weak impurities. As is done experimentally, 
we focus on the intensity in energy-momentum space 
of the Fourier-transformed differential tunneling conduc- 



tance, i.e., the magnitude of |j4^(q, J5)| in Eq. (21) 



1. QPI and renormalized band structure 

We start the discussion of QPI with a plot of 
(q,eV)| along = Qy for tf/tc = 0, Fig. Jgj to be 
compared to the single-particle spectrum shown m Fig. [6] 

To understand the QPI result, one has to recall that 
high QPI intensity at a given energy oo is expected - 
within a quasiparticle picture - for wavevectors which 
connect pieces of the quasiparticle iso-energy contour at 
UJ. Thus, for a simple inversion- symmetric dispersion ek, 
high QPI intensity at cj = ek will occur for q = 2k such 
that q connects the equal-energy momenta — k and k. 

Indeed, by using k = q/2 and "unfolding" the plot 
of |^J^(q, eV)| (i.e. adding a mirror image), one can 
qualitatively recover from Fig. |9] the single-particle spec- 
trum shown in Fig.[6j In particular, the reconstruction of 
the band structure upon cooling, where spectral weight 
is transferred within the optical gap near the Fermi level 
(Sec. VB ), is clearly visible in the QPI plot as well. Fig. [9] 

Of course, more complicated band structures, their 
nesting properties, and the effect of van-Hove singular- 
ities are not captured by the above argument. This is 
nicely visible in our data as well: The QPI signal near the 
middle of the band is influenced by the nesting properties 
of the square-lattice dispersion, Eq. (24). The iso-energy 



contour is a rotated square with perfect nesting in the 
non-interacting limit. As a result, that QPI response is 
very strong at q = (tt, tt), but extends to all wavevectors 
with Qx = Qy^ albeit with a reduced intensity, thanks to 
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FIG. 9: a) Intensity map of the Fourier transform of the spatially inhomogeneous part of the differential tunneling conductance, 
1^ J^(q, c<;)|, induced by a single point-like weak scalar impurity (i.e. the QPI signal), along the qx = Qy direction of the first 
Brillouin zone. The perculiar behavior near uj/W = 0.5 is due to perfect nesting of Ck, for details see text, b) Close-up near 
the Fermi level to be compared with the single-electron spectrum in Fig. [gJ). The model parameters are as in Fig. [3] the ratio 
of the tunneling amplitudes is tf/tc — 0. 



the interaction-induced quasiparticle broadening, com- 
pare Fig. [9|i near uj /W ^ 0.5. 



2. Temperature evolution of QPI 

An alternative way to represent QPI data and to 
discuss its relation to the single-particle spectra is via 
constant-energy cuts. Such plots are displayed in Figs.[l£ 
and [TT] which show the QPI signal together with the 
c and / spectra for different temperature and energies 

= and uo = A.7 x 10~^ (inside the hybridization gap), 
respectively. 

At the Fermi level, Fig. [Toj the c-electron spectrum il- 
lustrates the temperature evolution from the small Fermi 
surface at high temperature (Fig. [lOfDl) to the large 
Fermi surface at low temperature (Fig. [l0]35). The /- 
electron spectrum displays a sharp Fermi surface only at 
low T, where the intensity exceeds that in the c spectrum, 
emphasizing that the low-energy renormalized band has 
more / than c character. 

The QPI signal has a strong intensity at a wave- vector 
q essentially when the c spectrum exhibits a sharp quasi- 
particle peak at q/2. In particular, the data at high and 
low temperature. Figs. [lO^l and [10^5, reflect the cor- 
responding Fermi surfaces. However, a unique extrac- 



experiments, sub-atomic resolution allows one to obtain 
information beyond the first Brillouin zone, such that 
those ambiguities can be partially resolved. At interme- 
diate temperatures T ~ T^, the quasiparticle peak in the 
c spectrum dissolves, and consequently QPI response is 
very weak and diffuse, Figs. [10^2, a3. 

It is instructive to compare the signal at the Fermi 
level, Fig. [Toj to that at an energy inside the hybridiza- 
tion gap, Fig. [TT] Again, at high temperature, we have a 



sharp quasiparticle peak in the c spectrum, correspond- 
ing to the band forming a small Fermi surface, and, ac- 
cordingly, an intense QPI response at the wavevectors 
connecting portions of this iso-energy contour. Upon 
lowering the temperature, both the quasiparticle peak 
and the QPI intensity decrease and essentially disappear 
completely at the lowest temperature. This nicely reflects 
the absence of well-defined quasiparticles inside the hy- 
bridization gap, i.e., all intensity inside this pseudogap is 
incoherent. 



3. Tunneling paths and QPI 

To complete our survey, we show the QPI signal at 
the Fermi level, a; = 0, for a finite ratio of the tunneling 
amplitudes tf/tc in Fig. 12 As discussed for the d iffer- 



tion of the band structure from QPI may be difficult: ential conductance in the clean system (Sec. V C[ ), the 



Figs. 10 il and 10 i5 look rather similar, due to the fact 
that the small Fermi surface has an electron volume of 
Tic ~ 0.6 while the large Fermi surface has a hole vol- 
ume 2 — Tic — Uf ~ 0.56, i.e., both Fermi surfaces yield 
similar wavevectors for elastic scattering. In actual STM 



destructive (constructive) interference between the two 
tunneling channels for a positive (negative) ratio tf/tc is 
apparent through a reduction (enhancement) in the QPI 
signal at low temperatures, while it is mainly unchanged 
at high temperature. The effect of the destructive inter- 
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ference can be so strong that it essentially suppresses the 
QPI signal at low energies and temperatures, as can be 
seen for tf /tc = 0.6 in Fig. 
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VI. SUMMARY 

We have analyzed the temperature-dependent elec- 
tronic spectra in the periodic Anderson model, a paradig- 
matic model for heavy- fermion formation, using dynam- 
ical mean field theory (DMFT) with Wilson's numeri- 
cal renormalization group (NRG). Particular attention 
has been paid to the temperature evolution of the low- 
energy spectra and the crossover from a '"small" Fermi 
surface of conduction electrons to a " large" Fermi surface 
of composite heavy quasiparticles upon cooling. To make 
contact with STM experiments, we have further stud- 
ied the differential tunneling conductance, related to lo- 
cal spectral properties of the model, and its modulations 
arising from impurity scattering processes via quasipar- 
ticle interference (QPI). In particular, we go beyond the 
limitations of previous analytical approaches, based on 
a slave boson mean-field theory, by fully accounting for 
interaction-induced broadening of spectral features away 
from the Fermi level and incoherent scattering processes 
at elevated temperature. 

For the clean system, the local differential tunneling 
conductance is shown to display an asymmetric peak. 



similar to the one observed experimentally, whose shape 
and intensity, apart from material specific details, depend 
on the ratio of the tunneling amplitudes tc and tf into 
conduction-electron and /-electron states, respectively. 
(Previous theoretical studies found a hard gap near the 
Fermi level and recovered the experimentally observed 
peak only by an ad hoc addition of a phenomenological 
quasiparticle broadening.) For positive ratio t//tc, the 
effect of destructive interferences between the two tun- 
neling channels shows itself as a decrease of the differen- 
tial tunneling conductance intensity at low temperature, 
while it is mainly inefficient at high temperature. 

By studying at the momentum-dependent spectral 
functions at different temperatures, we unveiled the dy- 
namical reconstruction of the Fermi surface which hap- 
pens mainly within the optical gap. In this window of en- 
ergies, spectral weight at high frequencies is transferred 
from the bare conduction band to the edges of the renor- 
malized bands when the temperature is lowered down to 
zero. Islets of spectral weight at the Fermi level persist 
but become smaller and smaller with decreasing temper- 
atures until merging completely with the incipient band 
structure in the T ^ limit. This underlines the energy 
dependence of the crossover from a "small" to a "large" 
Fermi surface: this crossover, characterized by dissolving 
quasiparticle peaks, happens at a frequency-dependent 
temperature. 

Quasiparticle interference induced by impurity scat- 
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terers offers an opportunity to follow such a dynamical 
reconstruction, as one can at least partially reconstruct 
the band structure from the dispersive high-intensity fea- 
tures of the QPI response. However, we have also shown 
that this may be seriously hindered by the destructive 
effect of the interference between the two tunneling chan- 
nels (here for positive ratio tf /tc)^ which may lead to an 
almost complete extinction of QPI for certain energies. 

Our detailed study spans many issues that are of rele- 
vance for existing and forthcoming spectroscopic mea- 
surements on heavy- fermion materials. It would be 
particularly interesting to study the energy-dependent 
crossover between small and large Fermi surfaces advo- 



cated here. Layered heavy-fermion materials, e.g. of the 
so-called 115 family (CeCoIns and relatives) are most 
promising in this context, and corresponding experiments 
are underway.l^ 
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